Translationally invariant discrete kinks from one-dimensional maps 
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For most discretisations of the <j>'^ theory, the stationary kink can only be centered either on a 
lattice site or midway between two adjacent sites. We search for exceptional discretisations which 
allow stationary kinks to be centered anywhere between the sites. We show that this translational 
invariance of the kink implies the existence of an underlying one- dimensional map — F{<j)n). 

A simple algorithm based on this observation generates three different families of exceptional dis- 
cretisations. 

PACS numbers: 05.45.Yv, 63.20.Pw 



Since the early 1960s, the (fA-eqadXioTi, 
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has been one of the workhorses of statistical mechanics 
and quantum field theory 2]. Its kink solution, 
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together with the sine-Gordon kink, are the simplest ex- 
amples of topological solitons. More recently, interest has 
shifted towards the discrete 0^-theories 0,01, 



20„ + 4>n 



/l2 



+ /(0n-l>ri,0n+l), (3) 



and their solutions. Here h is the lattice spacing: 0„(i) = 
(j){x„,t), with Xn = hn, and the function / is chosen to 
reproduce the nonlinearity (^) in the continuum limit: 



(4) 



For a variety of discrete nonlinearities /, Eq.© admits 
stationary kink solutions 0, 0. The discrete (/)^ kinks 
have been used to describe incommensurate systems and 
narrow domain walls in ferroelectrics and ferromagnets, 
topological excitations in biological macromolecules and 
hydrogen-bonded chains, and polymerization mismatches 
in polymers 0, Q • Physically, one of the most signifi- 
cant properties of domain walls and topological defects is 
their mobility however whether the discrete equations 

fcan admit travelling kinks remains an open question 
II B E Ei El- The continuous (/i^-equation Q is 
Lorentz- invariant, and so the existence of the travelling 
kink ^ is an immediate consequence of the existence of 
the stationary soliton, Eq.(|2Jl with c — 0. The discreti- 
sation breaks the Lorentz invariance and the existence of 
travelling discrete kinks becomes a nontrivial matter. 

In fact, the discretisation even breaks the translation 
invariance of Eq.l^. As a result, the stationary kink 



can be centered only at a countable number of points — 
usually on a site and midway between two adjacent sites 
1^ Q . This breaking of the translation invariance is con- 
nected with the presence of the Peierls-Nabarro barrier, 
an additional periodic potential induced by discreteness. 

Miraculously, there are several exceptional discretisa- 
tions which, while breaking the translation invariance of 
the equation, allow the existence of translationally in- 
variant kinks; that is, kinks centred at an arbitrary point 
between the sites. One such discretisation was discovered 
by Speight and Ward using a Bogomolny-type energy- 
minimality argument p, 10]: 
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Another one derives from the Ablowitz-Ladik integrable 
discretisation of the nonlinear Schrodinger equation; it 
was reobtained by Bender and Tovbis [13j from the re- 
quirement of suppression of the kinks' resonant radiation: 



,/ - 4 {^n+l + K-l) (1 - eg) . 

Finally, the nonlinearity 
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was identified by Kevrekidis 12), who demonstrated the 
existence of a two-point invariant associated with the sta- 
tionary equation 



(f)n+l - ^(fln + (f>n- 



+ /((/.„_!, 0„, )= 0, (8) 



with / as in jnj and 0. 

Although the translation invariance of a stationary 
kink docs not automatically guarantee the existence of 
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a travelling soliton, it is natural to expect it to be a 
prerequisite for kink mobility. For example, in the vari- 
ational description of the slowly moving kink, the solu- 
tion is sought as a stationary kink with a free continuous 
parameter defining its position on the line |^. Also, the 
Stokes constants measuring the intensity of resonant radi- 
ation from the translationally invariant kinks were found 
to be at least an order of magnitude smaller than the cor- 
responding constants in models with noninvariant kinks 
With an eye to a future attack on travelling kinks, 
it would be useful to identify all discretisations of the 
(j)* theory supporting translationally invariant stationary 
kinks. The purpose of this note is to provide a general 
recipe for the generation of such exceptional discretisa- 
tions /(0„_i, (?!)„, 0„+i). 

We start with a simple observation which, however, 
holds the key to our construction. Assume we have a 
nonlinearity / which supports a stationary discrete kink, 
i.e. a monotonically growing sequence (/)„: — 1 < 0„ < 
(j^n+i < 1 for — cxD < n < oo. Furthermore, assume there 
exists a continuous monotonically growing function g{x), 
defined for all real x, such that = g(n). (The func- 
tion g can also depend on h parametrically but we omit 
this dependence for simplicity of notation.) This func- 
tion generates a family of kinks centered at an arbitrary 
point x'"^ on the line: = g{n — x^^">). It is important 
to emphasise that such a continuous function can exist 
only for exceptional discretisations. For generic discreti- 
sations, the function g{x) can be defined only on integers. 

The existence of the function g(x) defined on the en- 
tire real line — or, equivalently, the translation invari- 
ance of the kink — implies that the stationary equation 
derives from a two-point map. Indeed, since g{x) is 
monotonic, we can write n = g^^{(f)n). Now since 4>n+i = 
.g(n-t-l), wehave(/)„+i = 5(3-1 -|- 1) = F(0„)^ which 
is a well-defined one-dimensional map. 

This observation suggests the following strategy for the 
construction of exceptional discretisations. Assume we 
have a ID map which we will write in the form 

- 4>n^ hH{(j)n+l,(j)n)- (9) 

Let H satisfy the following continuity condition: 

i/(0,0) = 1(1-02), (10) 

This condition is necessary to make sure that the map 
becomes 

= (11) 

in the continuum limit. The stationary (c = 0) kink 
solution (O of Eq.(^ is, simultaneously, a solution of the 
first-order equation (|ll|l . Imposing l|l()(l we ensure that 
the discrete kink of Q will have the correct continuum 
limit. Next, Ec ljlUII implies that the map © has just 



one pair of fixed points, 0* = =Fl- For small h, (j)n+i 
remains close to 4>n and hence, i?((/)„+i, 4>n) remains close 
to l(TT)|l which is positive for < 1. Consequently, no 
matter what j^ol < 1 we start with, the sequence (j)n is 
monotonically growing — at least until is not very 
close to 1. To ensure that it remains such near the fixed 
points, we assume that 0* = — 1 is a source and 0* = 1 a 
sink. (That is, small perturbations S(f>n = 0„ — 0* satisfy 
S(t>n+i = ^S(f)n with A > 1 near (f>^ = —1 and < A < 1 
near (j)^ — 1.) Then, for any h smaller than some h and 
any (j)Q between —1 and 1, there is a number TV such that 
\(t>N — 4'* \ is so small that all (f>n with n > N are entrapped 
by the "linear neighbourhood" of 0* = 1 and those with 
n < —N are all in a neighbourhood of = —1. This 
means that each (f>Q with |(/)q| < 1 defines a monotonic 
kink solution and so for any sufficiently small h we have 
a one-parameter family of stationary kinks. Speight \vi\ 
gives a less intuitively appealing but more rigorous proof 
of this fact; he also shows that our assumption on the 
character of the fixed points can be relaxed. 

Next, squaring both sides of © and subtracting the 
square of its back-iterated copy, 

4>n - 4>n-l = hH{(j)n, (j)n-l), (12) 

produces an exceptional stationary Klein-Gordon equa- 
tion 

h'^ Ipn+l - 4'n-l 

(13) 

If H is symmetric: H{(f)n, 4>n~i) ~ H{(f)n-i, 4'n), the nu- 
merator vanishes exactly where the denominator equals 
zero, so the discretisation ()ll-t|l is nonsingular. 

If we want to have polynomial discretisations of the 
(j)* theory, the function has to be a quartic polyno- 
mial. This leads to two possibilities, one where H is the 
square root of a polynomial, and the other where if is a 
polynomial itself. These can be written jointly as 

(0„+l - 0„)" = /i'"P2m (</-«+!>«), (14) 

where m = 1 or 2, and P2m{u, v) is a polynomial of degree 
2m that satisfies the symmetry and continuity conditions 

P2m{u,v) = Plm{v,u), (15) 
P2raM) = 2-™(l-02)". (16) 

The condition H16() is a consequence of Ea. (|10|l . 

Before we proceed to the classification of the resulting 
models, it is pertinent to note that the linear part of the 
function / in iPJ can always be fixed to without loss 
of generality. Indeed, the most general function satisfying 
Q is / = a0„ + i (i - a) (^n-i-i + 0n-i) + cubic terms. 
Since in (|S1) is a free parameter, we can always make 
a replacement h h such that a — = d — 2/h'^. In 
particular, we can set a = ^ which gives 

f{4>n-lAnAn+l) = T^'Pn - Q{4>n-1, (l>n, 4'n+l) , (17) 
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where Q is a homogeneous polynomial of degree 3. 

Let, now, m = 2 in Ea. (|14(l . Provided P4 satisfies 
conditions l(T^ and ((TCI) , the numerator -^^((^n+i, ~ 
PA{4'ni4>n-i) of the fraction in the right-hand side of 
Ea. ((TC|l divides ((/)ti+i — (t^n-i) and so Ea. ((TC|) will be 
of the form © with some cubic function /. The most 
general choice for such a polynomial is 

Pi{u, w) = i - - vf - ]^uv 

+ ^ + '"^) + Puv{u^ + v^) + ju^v^] , (18) 

where a, (3, 7 satisfy 2a + 2j3 + ^ — b and ^ is arbitrary. 
Picking the positive value of ^/Pl and assuming that h 
is sufficiently small, one can check that the fixed points 
= =pl of the map l|14|) are a source and a sink, for any 
/i, a and /3. Consequently, the resulting cubic polynomial, 

<3 = ^ [a{(t)n+l + (l)n-l){4>'i+i + (fi-i) + 70^(0n+l 
+ 0„_i) +(3(j)n{4>l+l + 0' + + 0«+10„-l)] (19) 

with 7 = 5 — 2(a + /3) defines a two-parameter family of 
models with translationally invariant kink solutions. 

The discretisation (|19|l includes, as particular cases, 
the Bender- Tovbis function © (which results from set- 
ting a — 13 = Q) and the Kevrekidis nonlinearity iQ (for 
which a — |, (3 — 0). Another simple function arises by 
letting a = 7 = 0; this is a new model: 

Q = ^0n(0^+i +<fin+ <i>n^l + 0„+l(?!)„-l)- 

Now let m = 1. The most general quadratic P2 satis- 
fying (dSJ-CSI) is 

P2{u,v)^'^-a{u-vf -'^uv, (20) 

with an arbitrary a. Note that for h < 2 and any a, — 
=pl are a source and a sink. Hence all resulting models 
will exhibit continuous families of kinks. Substituting 
Ea. (|20|l for H in (|13(l . we obtain just a particular case 
of the nonlinearity (|19|l . corresponding to the choice of 
the quartic (|18|l in the form of a complete square: P4 — 
P2 ■ To obtain new models, we need to note a symmetry 
/3(^„_i, (/)„, 0„+i) = which follows from Eq.lQ^ with 
m ~ I. Here 

-^3 = P2i4>7i+l,4>n){4'n~(t>n^l)+P2i4'n,4'n-l)i4'n-<t'n+l)- 

Equation H13|) remains valid if pi^ is added to its right- 
hand side, with an arbitrary coefficient /3. The resulting 
function Q has the form 

Q = a^(</'n+l + </'n-l) + 27 (a - (3) Ipn+lMn-l 

+ a{a - 13) (?!)„+i0„_i(0„+i -f- 0„_i) 
+ [2a2 -f 72 + /?(7 - a)] (pliK+i + K-i) 
+ a (27 + /3) 0„(02^i + <Pl_,) + 2a (7 + /?) </>L (21) 



where 1 = \ — 2a. Ea. (|21|l defines a two-parameter fam- 
ily of discretisations supporting translationally-invariant 
kinks. These models cannot be obtained within the 
m = 2 analysis above — unless /3 = 0, of course. 

Letting a — (3 = ^, we recover the model of Speight 
and Ward, Eq.(|SJ|. Another particularly simple, new, 
model is obtained by taking a = and (3 = 

It is instructive to compare discretisations furnished 
by our method with those arising from the requirement 
of the absence of resonant radiation from the kink [isj . 
The advance-delay equation associated with Eq.©, 

(j){x + h)- 2(j){x) + (l){x ~ h) 

+ h^f{(t){x-h),(l){x),<j>{x + h)) = 0, (22) 

can be solved to all orders as a perturbation expansion in 
h |l4| ; the resultant solution depends continuously on the 
position parameter x^^\ It is, therefore, only the terms 
which lie beyond all orders in h which present an obsta- 
cle to the existence of a translationally invariant kink. 
These terms vanish if Stokes constants at all orders van- 
ish. The leading Stokes constant vanishes if there exists 
a convergent solution in powers of to the equation 

ip{z + 1) - 2ip{z)+ip{z - 1) 

- Q{ip{z - l),ip{z),ipiz + 1))|^^Q = (23) 

(see e.g. ^3)- This equation comes from a rescaling of 
Eq. (|22|l near the singularities of its leading-order solution 
^ at Xn = TTi{l + 2n), n E Z, in the limit h 0. 
The convergence of a power-series solution to Eq. (|23|l is 
necessary for the absence of oscillatory radiation tails in 
its "parent" equation (|22|l . 

In general, a numerical procedure is required to deter- 
mine whether the series converges for a given Q, but we 
can easily generate a class of models for which it trun- 
cates after the first term. This was the method employed 
in Ref. 13] in deriving Eq.ljHJ. It is a matter of direct sub- 
stitution to check that the most general cubic polynomial 
for which <y9 = 2/z is a solution of Eq. (|23|) , is 

Q = Cr(j)n{(t>n+1 + " 2a(j)n+l(t)n-l{4'n+l + (?!)„-l) 

+ Q - fpliK+l + (l)n-l) + /3(/i„_i0„ (/)„+!, (24) 

with a, (3 arbitrary constants. 

The fact that the leading-order Stokes constant is zero 
is necessary but not sufficient for equation (|22|l to have 
continuous families of kinks for finite h. We tested, nu- 
merically, a particular representative from the class H24|l : 

Q = <Pn+l(Pn-l 2 4 ■ ^^"^f 
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This model is obtained by letting a = —j and /? = i. 
To check whether translationally-invariant kinks exist 
or not, we have computed a stationary on-site kink for 
the model (|25|l and calculated eigenvalues of the associ- 
ated linearised operator for an equidistant sequence of 
/i- values, ranging from h = to h = 1.179 with in- 
crement 0.001. For h smaller than 0.556, the smallest- 
modulus eigenvalue was found to be smaller than 10~^^ 
which is our numerical error of computation. However 
as h increases from 0.556, the smallest eigenvalue grows 
to A = 9 X 10~^ aX h = 0.955, then decreases, crosses 
through zero at h = 0.993, after which grows in modu- 
lus to A = -6 X 10-4 at h = 1.179 (FigHJ. Thus, the 
zero eigenvalue, indicating the existence of a continuous 
family of solutions, is not present in the spectrum for 
h > 0.556. For h < 0.556, the smallest A is apparently 
also nonzero (though very small). The only exception is 
the value h — 0.993 for which the zero mode does exist. 
This means that the model H25(l supports a continuous 
family of kinks for just one, isolated, value of h. 

However, there exists a family of exceptional discreti- 
sations which reduces to (|24|l in the limit h ^ 0. Indeed, 
when a = in Ea. (|20() . the map (|14|l has one more sym- 
metry: /2(0„_i, 0„, 0„+i) = 0, where 




We can add a{(j)n+i + (j>n-i)l2 to the right-hand side of 
(tT!^ . along with f]!^ (with cr, (3 arbitrary constants.) This 
gives rise to the following family of discretisations: 

Q = Cr<j)ni(l>n+1 + 4>n-lf - 2a-</)„-10„+l ((/)„+l -|- 4>n-l) 




+ -Crh^(f,n{(l)n+l+(l)n-lf- (26) 

Except for the last, 0{h^), term, this coincides with 
Eq.(j2i|). 

For the map (fnjl + (pn|l with a — 0, the kink can be 
found explicitly. This implies that the discretisations 
(I26|l also share an explicit kink solution (for all f3 and 
a): (pn — tanh(art — x^"-*), with tanha = h/2. 

Our final remark is on the conserved quantities of 
Eq. (PJ . The translation invariance of the stationary kink 
does not imply the invariance of equation Q and hence 
the conservation of momentum. The discretisation H13|) 
[and hence H19() ] conserves momentum whereas the 
nonlinearities H21|l and (|26|l (with /3,cr 7^ 0) — do not. 
Moreover, that the discretisation / is exceptional does 
not guarantee that Eq.Q has any integral of motion 
whatsoever. In particular, out of the three families (|19|1 . 
(|21|l . and H26() . only one model conserves energy, namely 
Speight and Ward's, Eq.©. 



In conclusion, we have identified three families of dis- 
cretisations of the 0^ equation which support translation- 
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FIG. 1: The smallest-modulus eigenvalue as a function of h. 
The cusp occurs at the point h = 0.993 where A changes sign. 



ally invariant stationary kinks: Eas. (|19|l . (|21|l and l|26(l . 

In each case we have exhibited, explicitly, the underlying 
ID map. 

LB. is a Harry Oppenheimer Fellow. 0.0. was sup- 
ported by a Nelli Brown Spilhaus scholarship. 



[1] A. R. Bishop and T. Schneider (eds). Solitons and con- 
densed matter physics. Springer- Verlag, Berlin, 1978 

[2] R. Rajaraman, Solitons and Instantons. North-Holland, 
Amsterdam, 1982 

[3] Yu.S. Kivshar, A. Sanchez, L. Vazquez, Phys. Rev. A 45, 
1207 (1992) 

[4] S. Flach, CR. Willis, Phys. Rev. E 47, 4447 (1993); S. 
Flach, K. Kladko, Phys. Rev. E 54, 2912 (1996); P. Ma- 
niadis, G.P. Tsironis, A.R. Bishop, and A.V. Zolotaryuk, 
Phys. Rev. E 60, 7618 (1999); P.G. Kevrekidis, M.I. We- 
instein, Physica D 142, 113 (2000); A.B. Adib, C.A.S. 
Almeida, Phys. Rev. E 64, 037701 (2001) 
[5] P. Prelovsek and I. Sega, J. Phys. C 14, 5609 (1981) 
[6] J. A. Combs and S. Yip, Phys. Rev. B 28, 6873 (1983) 
[7] P. Bak, Rep. Prog. Phys. 45, 587 (1982); M. 
Remoissenet, Waves Called Solitons. Springer- Verlag, 
Berlin, 1996; O.M. Braun, Yu.S. Kivshar, Phys. Rep. 
306, 1 (1998) 

[8] S. Flach, Y. Zolotaryuk, and K. Kladko, Phys. Rev. E 
59, 6105 (1999) 

[9] J.M. Speight and R.S. Ward, Nonlinearity 7, 475 (1994) 
[10] J.M. Speight, Nonlinearity 10, 1615 (1997) 
[11] J.M. Speight, Nonlinearity 12, 1373 (1999) 
[12] P.G. Kevrekidis, Physica D 183, 68 (2003). 
[13] CM. Bender, A. Tovbis, J. Math. Phys. 38, 3700 (1997) 
[14] O.F. Oxtoby, D.E. Pelinovsky, and I.V. Barashenkov, 

submitted for publication 
[15] A. Tovbis, M. Tsuchiya, C. Jaffe, Chaos 8, 665 (1998) 



